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We study the long-range Coulomb interaction effects on the double-Weyl fermion system which is possibly 
realized in the three dimensional semimetal HgCr 2 Se 4 in the ferromagnetic phase. Within the one-loop renor¬ 
malization group analysis, we find that there exists a stable fixed point at which the Coulomb interaction is 
screened anisotropically. At the stable fixed point, the renormalized Coulomb interaction induces nontrivial log¬ 
arithmic corrections to the physical quantities such as specific heat, compressibility, the electrical conductivity, 
and the diamagnetic susceptibility that are obtained utilizing the renormalization group equations.. 


I. INTRODUCTION 

There has been recently much interest in semimetals, which 
support gapless quasiparticle excitations only in the vicinity 
of isolated band touching points in the Brillouin zone (BZ). 
When the Fermi energy is pinned to the band touching points, 
these semimetals can possess universal power-law behaviors 
for thermodynamic and transport quantities as a function of 
temperature or external frequency. There are many well 
known experimental examples of the semimetals which pos¬ 
sess linearly dispersing massless Dirac quasiparticles in both 
two dimensions ( 2 D) and three dimensions ( 3 D), e.g. Mono- 
layer graphenei— in 2 D and Bii-ajSbj?^— , Pbi_a;Sn2;Te^, 
and Cd3As2^, NaaBiiS in 3 D. It is also possible to realize 
parabolic semimetals which possess parabolic dispersions at 
band touching, e.g. Berner-Stacked bilayer grapheneii in 2 D, 
and HgTei, gray tii>i^, and the normal state at high tempera¬ 
ture for some 227 irradiates such as Pr2lr207ii^^ in 3 D. 

In the presence of strong spin-orbit interactions in three di¬ 
mensions, the unusual semimetallic phase called the topolog¬ 
ical Weyl semimetals may exist and have been confirmed in 
TaAs recentlyi^— . The Weyl semimetals are also predicted to 
exist in pyrochlore iridates^^^, cold atom systems2L2Z, and 
multilayer topological insulator systems22i22. In close prox¬ 
imity of the gapless points, the effective Hamiltonian is de¬ 
scribed by a two-component wave-function termed the Weyl 
fermion and the gap closing point is the Weyl node. The Weyl 
nodes are protected from opening a gap against infinitesimal 
translations of the Hamiltonian; these points act as monopoles 
(vortices) of 3 D Berry curvature, as any closed 2 D surface sur¬ 
rounding one of them exhibits a finite Chern flux, and a Weyl 
node can only be gapped by annihilation with an anti-Weyl 
node of opposite monopole charge. 

In addition to the usual Weyl semimetals, recently Ref. 
proposed the possible presence of new 3 D topological 
semimetals in materials with point group symmetries termed 
as double-Weyl semimetals. The new double-Weyl semimet¬ 
als possess Weyl nodes with quadratic dispersions in two 
directions, e.g. x — y plane. The double-Weyl nodes are 
protected by C4 or Cq rotation symmetry and are suggested 
to be realized in the 3 D semimetal HgCr2Se4 in the ferro¬ 
magnetic phase, which possess a pair of double-Weyl nodes 
along r Z direction^SiiL. The first-principle calculations in the 
material HgCr2Se4^ also suggested the existence of double- 
Weyl nodes, which is qualitatively in agreement with the re¬ 


cent transport experiments in HgCr2Se4^ that confirm the 
half-metallic property of the HgCr2Se4. The (anti-)double- 
Weyl node possesses a monopole charge of (- 2 )h -2 and the 
double-Weyl semimetal shows double Fermi-arcs on the sur¬ 
face BZ^Sii. This new semimetallic phase with an in-plane 
quadratic dispersion can serve as a new platform for studying 
the (long-range) Coulomb interaction effects on the double- 
Weyl fermion. The low-energy physics of the double-Weyl 
fermion can possibly serve as a new source term contributing 
to the physical properties in HgCr2Se4 for chemical potential 
sitting near the Weyl nodes, such as a term to the specific 
heat C that was not considered previously^. 

In this work, we consider a single double-Weyl fermion 
coupled to the long-range Coulomb interaction and study the 
effects within the one-loop renormalization group (RG) anal¬ 
ysis in the Wilsonian momentum shell scheme^. Due to the 
anisotropic dispersions of the double-Weyl node, the den¬ 
sity of states (DOS) is linearly proportional to the energy, 
D{e) oc e, sharply different to that in the usual Weyl fermion 
with D{e) oc e^. Due to the anisotropic dispersions, the scal¬ 
ings for the three spatial coordinates can be different. In the 
noninteracting limit, for scaling transformation invariance of 
the action we And the dynamical scaling exponent z = 2 , the 
scaling exponents of the spatial coordinates x and y are 1, i.e. 
the scaling dimension [x] = [^ = — 1, while that of Zis two, 
[2] = — 2 . The result of such nontrivial scaling transforma¬ 
tions in spatial dimensions result in the Coulomb interaction 

with engineering scaling dimension [e^] =2 — 1 = 1 and 
anisotropy parameter p, which dictates the anisotropy of the 
system, with enginnering scaling dimension [p] = — 2. 

After coarse-graining within RG analysis, we find that in 
the low-energy limit, the system becomes highly anisotropic 
and 77 becomes infinitesimal. The Coulomb potential receives 
strong renormalization along the linear dispersion axis and 
the Coulomb interaction strength also becomes infinitesi¬ 
mal due to the strongly irrelevant anisotropy variable p. The 
composite variable similarly equal to the ratio of Coulomb 
interaction strength and the square root of the anisotropy pa¬ 
rameter, ^ e'^/y/rj, approaches & fixed value in the low-energy 
limit, which deflnes the stable fixed point in RG. At the stable 
fixed point, the Coulomb potential is renormalized anisotrop¬ 
ically, which is consistent with the random phase approxi¬ 
mation (RPA) calculation detailed in the supplemetal mate¬ 
rial. Furthermore, we And that the square of the Coulomb 
interaction,~ e^, under coarse-graining process decreases in 
a logarithmic manner, which leads to logarithmic suppres- 
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sions to several physical quantities such as specific heat C 
compressibility k and the so called finite frequency (dynamic) 
conductivity <J^^=x,y,z(!^), and the dc conductivity 
Furthermore, we find unexpectedly that the diamagnetic sus¬ 
ceptibility xd gets enhanced in a logarithmic manner due to 
the long-range Coulomb interaction. 

The paper is organized as follows. In Sec.|II]we introduce 
the model Hamiltonian followed by weak-coupling RG analy¬ 
sis. In Sec.|III]we utilize the RG equations near the fixed point 
to obtain logarithmic corrections to various physical quatities. 
In Sec.|IV]we conclude with some discussions. 


II. DOUBLE WEYL SEMIMETAL IN THE LONG-RANGE 
COULOMB INTERACTION 



We consider a minimal model of a single double-Weyl 
fermion coupled to the long-range Coulomb interaction. The 
action in the Euclidean path integral formalism is 


Sl= drd^xl '0^ dr — iecj) + d(—iV) • a ip + 


[{dxCpf + {dy(pf] + ^( 5 ^ 0 )^|,( 1 ) 


with d = {-m~^ ipl - dy) , - m~^{ 2 dxdy), - iv^dz}, 

where ip and (p represent the electron annihilation field and the 
boson annihilation field. The integration of (p gives the usual 
instantaneous long-range Coulomb interaction. The variable 
m along the x — y plane represents the effective mass and the 
variable Vz along z direction is the velocity. The anisotropy 
variable p is introduced due to the anisotropic dispersions. 

We choose the scaling transformations for the fields and 
the variables as r = x = b Xr, y = b y^, z = 

b""^ Zr, Vz = Z-^Vz,R, m-^ = p = Z-^yR, 

e = Ze^^'^CR, Ip = ZpjP^^'^ipR, 0 = ZpjP^^'^cpR, where the pa¬ 
rameter b = represents a length scale slightly greater than 
one with a logarithmic length scale £ <C 1 and the subscript R 
labels renormalized variables during coarse-graining. 

For analytically extracting the RG equations, we adopt 
the RG scheme in Ref. [HI and perform integration within 

a momentum shell in the q± = ^ direction, i.e. 

q±^ € [Ae“^, A], and no restriction along qz direction {\qz\ € 
[0, Az —>■ oo]). For clarity in the presentation of RG results, 
we introduce two variables a = and A = 

We obtain the one-loop RG equations, detailed in App.|A] 
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FIG. 1. RG stream plot for the model action, Eq. ([TJ. The red dot 
represent the fixed point (a*, A*) = (0,1/2), where the Coulomh 
interaction receives anisotropic screening. The RG streams mostly 
flow to the stable fixed point along the A = 1/2 path. 


where c = ln(3 -F 2y/2)/2 ~ 0.881. If we hold Vz and m ^ 
fixed, we get 


z = 2 - 


Zi = z + 


3c a^ 

TT’ 

3^2 

iT 


= 2-F 


3(1 — c) 

4 T 


( 6 ) 

(7) 


We can see that fixed points are located at (a, A) = (0, 0), 
and (0, 1/2). Linearizing around these two fixed points, we 
find that the fixed point (0, 0) is the unstable Gaussian fixed 
point, and (0, 1/2) is the stable fixed point controlled by the 
parameters a and A. The RG streamplot is shown in Fig. [T] 
The Coulomb interaction decreases to the stable fixed point 
mostly along the path of A = 1/2. Along this path, the 
square of the Coulomb interaction (e^) decreases to zero in 
a logarithmic manner, which is reflected as the nonmonotonic 
temperature or frequency dependences of physical quantities 
such as the specific heat C, compressibility k, the finite fre¬ 
quency (dynamic) conductivity a'y.y.=x,y,z{'^), the dc conduc¬ 
tivity and the diamagnetic susceptibility xd, which we 
will show below. 

At the fixed point (cts, A^) = (0,1/2), the dynamical ex¬ 
ponent z = 2 = zi- If we focus on the renormalized boson 
propagator (which gives the Coulomb screening effects), at 
the fixed point it is (below we will suppress the irrelevant di¬ 
mensionful variables) 


9j. + 92 ~ ^^(q) ^ 9l(l + ^s(-) + 92(1 + 4As£) 

= ql+ql{l + 2£), (8) 


which shows that there is only a correction along the qz direc¬ 
tion. If we consider the relative scalings between q± and qz, 
we can obtain 

9! + 92 ~ 11(9) 9i. + 9219^1 =9i + l92|, (9) 


where zi = 2 at the stable fixed point. The renormalized 
Coulomb interaction at the fixed point becomes V/(q) ^ 
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\/{q\ + |9z|)- Fourier transforming Vc{q) back to the real 
space, we find that the renormalized Coulomb potential falls 
off anisotropically, Vc{r^,\z\ = 0) ^ and Vc(r± = 
0, l^l) ^ where r± = + y^. In the Add. IbI we 

perform RPA analysis of the screened Coulomb interaction 
and we find that the scaling analysis above is consistent with 
the RPA analysis. 


III. LOGARITHMIC CORRECTIONS TO THE SCALING 
BEHAVIOR OE PHYSICAL QUANTITIES 


According to the strong-coupling analysis in App. we 
find that the renormalized Coulomb interaction causes an in¬ 
frared logarithmic divergence that can lead to the logarith¬ 
mic corrections to physical quantities. Instead of calculating 
higher-order corrections in the perturbation theory, we follow 
Ref. to utilize the RG equations near the stable fixed point 
and the scaling arguments to obtain the scaling behaviors of 
the physical quantities. 

Focusing on the path of A = 1/2 near the stable fixed point, 
we know that the mass inverse m~^ and velocity receive 
nontrivial renormalization as m~^{l + In b) and Vz{^ + 
In6), where \nh = ^. From the scaling invariance of ac¬ 
tion, we know Z^-i = (l + In 6), 

and Zy^ = (l -f In 6). For a finite potential term, 

we can consider to add a term —p / (fix f to the ac¬ 

tion and we can obtain the transformation jj, = b ~^We 
can also consider the free energy F, which can be a general 
function of T, a, p, magnetic field B, etc., that transforms 
as F = Zp^Fr. Considering the exponent of a partition 
function J fix J drF, we know that the exponent should be 
dimensionless, which requires F = Fr. In ad¬ 

dition, the temperature transforms under coarse-graining as 
T = b-^TR. 

The relevant RG equations near the stable fixed point along 
the path A = 1 /2 are 


da 
din 6 
dT 
din 6 


3c-f 2 


a^(6), 


T{b) 



( 10 ) 

( 11 ) 


Solving the RG equations, we get 

a^(6) = (l-I-(3c-|-2)a^ln6) (12) 

T{b) = Tfi (1 -f (3c-f 2)a2 In 

ya^lnfo^ , (13) 


The renormalization of specific heat under RG can be obtained 
as C = —'rfjS = Choosing the cut-off b to 

be b* and using the high temperature result of specific heat, 
Cr ^ Tq which is originated from D{e) ~ e in the noninter¬ 
acting double-Weyl semimetals, we get 

rp2 rp2 

(1-f |a2 In (1-f f a2 In + 3^2 > 

(15) 

where we crudely approximate c ~ 1 in the last line. The 
compressibility can be obtained as k = kr. 

If we substitute b* for b and use the noninteracting result kr ~ 
To, we get 


T 

(l -I- In -I- In 


T 

(l -I- In^)^ 

(16) 


Let’s shift our focus on the scaling behavior of the fi¬ 
nite frequency conductivity and dc conductivity. In the 
presence of magnetic vector potential, the kinetic terms are 
modified as —idj —>■ —idj + eAj, with j = x, y, z. 
Due to the minimal coupling, we require that the compos¬ 
ite variables eAj rescale the same as that of dj, which leads 

to A_l = 6“^ [l-f In 6] [l-f 2 q;^ ln6] A^^r, 
with A± = A^/y, and Az = 6“^ [l-I-ln6] x 

X [l -f |q;^ In 6] ^ [l + 2 q;^ In 6] Az^r, where we ex¬ 
plicitly use A = 1/2 near the fixed point and the scaling of 
electric charge Ze obtained in the RG equations derivation in 
App.m 

In order to obtain the scaling relations for the conduc¬ 
tivity, we can rely on the Kubo formula. According to 
the Kubo formula, the current-current correlation function 
App{iuJn) = f dre**^”^ / d^g (T^ [jy(q, 0)]) 

can be related to the dynamic conductivity as 
<Tfj.fj.{uj) = -ImIly,p{oj)/u}, where w -f i0+. j±{q,T) 

and jz ((/, t) can be obtained from the Fourier transform of 
j±ix,T) ~ em-'^'ip^x,T){-iV) ■ r) and/^(f,r) = 

evz^p"^{ x,t)(7^'iI){x^t). We first obtain that j±{q,T) = 

b~^ [l + \nb] [l -f 2 a‘^ In b] 2 x,r and jz (t t) = 
[l -I- In [l -I- 2q;^ Ind] which lead to 

aXX = 6“^ [l + Infe] ^ [l-f 2a^ In 6] crand 

Ozz = [l + In 5] ^ [l + 2 q;^ In b] azz,R- For uj > T 
(with T —0, w —?► 0), we introduce the cut-off frequency 
uiQ = x!{b*) with 


where a and T represent the bare Coulomb interaction a(0) 
and bare temperature T(0). Choosing the temperature cut-off 
T(6*) = To = we get 



3 c 2 1 ^0 

1 -I- —a In — 
4 u) 


(17) 


b* 


n 

T 



3 c 2 n Tq 
—a In — 
4 T 


which is due to the fact that the uj{b) scales the same as the 
temperature T(6), i.e. O = b~^OR, where O = T, oj. With 
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the cut-off frequency, we obtain 




[1 + T«^ln(i^)] [l+3„2in(^^)] 

3c -I- 1 2 T 
1- :—a In 




o-zz(w) 




e) 




e^uj, 


1-In 

4 


B) 


(18) 


e^. 


(19) 


where we approximate crj_j_^/j(a;) and (T^_h(u;) to be the non¬ 
interacting results of the dynamic conductivity at finite fre¬ 
quency obtained in App. |D] 

The dynamic conductivity calculations at p <C T in the 
noninteracting limit in App. |Dj also show the existence of 
the Drude peak and the linear-T dependent xx /j/y-component 
dc conductivity, and T-independent zz- 

component dc conductivity, Following similar discus¬ 
sions above with high temperature cut-off, Eq. (O, the dc 
conductivity also receive nonmonotonic temperature suppres¬ 
sion, which are similar to Eqs (fT8l) -(fT^ with uj ^ T, 


ai\{T) 




[1 

^ToY 

[ 4 



e^T, 


l--a^ln(^ 


( 20 ) 

( 21 ) 


Last but not least, we focus on the temperature dependence 
of the diamagnetic susceptibility. The diamagnetic suscep¬ 
tibility can be obtained from taking second derivative of the 
free energy, xd = —d^F/dB^. Since B = V x A, the renor¬ 
malization of the magnetic field under coarse-graining can 
be obtained straightforward using the renormalization of Aj 
illustrated above. We find that the diamagnetic susceptibility 
scales differently for in-plane magnetic field B = Br± and 
for perpendicular magnetic field B = Bz. Eor B = Bf±, 
the diamagnetic susceptibility renormalizes as Xd — 

[l -I- In b] [l -I- \c? In b] [l -|- 2 q;^ In b] Xd,r- 
Eor B = Bz, the diamagnetic sus¬ 
ceptibility renormalizes as xh = 

6“^ [l -I- In b]^ [l -f |q;^ In b] ^ [l + In b] xh r- 
We use Eq. (fT^ and the noninteracting results of xd,r 
which we derive using the Eukuyama formula for the orbital 
diamagnetic susceptibility^ in App. |E] We obtain 




3c -f 5 2 


^To 

Xd ^ 

1 -f 

4 “ 

ln| 

kT 



6c-2 2 


^To 

Xd - 

1 + 

4 “ 

ln| 

kT 




e^T 


mvz 


( 22 ) 

(23) 


We expect that the temperature dependence of the diamagnetic 
susceptibility for a magnetic field in general direction at low 
temperature to be 


XD 


1 -f In 



(sin^ 0 xo + cos^ 0 T) , (24) 


where 0 is the angle between the magnetic field and z-axis, 
i.e. B ■ z = B cos 9, and xo is ^ constant independent of 
temperature and the diamagnetic susceptibility gets enhanced. 


IV. DISCUSSIONS 


We study the long-range Coulomb interaction effects on 
the double-Weyl semimetals. Within one-loop renormaliza¬ 
tion group analysis we find that the composite variable de¬ 
fined as the ratio of the Coulomb interaction strength and the 
square root of the anisotropy parameter, ~ e^/y^, is fixed 
to be finite at long-wavelength, which defines the fixed point. 
Eocusing near the fixed point, we utilize RG equations to ob¬ 
tain nonmonotonic temperature or frequency dependences of 
various physical quantities. 

Though the long-range Coulomb interaction induces log¬ 
arithmic corrections to several physical quatities in experi¬ 
ment, the fundamental Berry curvature structure around the 
double WeyTpoint remains unaltered, similar to the situa¬ 
tions in single-Weyl point and the Dirac points of graphen^. 
In the presence of Coulomb interaction, the renormalized 
low-energy description near a double-Weyl point is similarly 
Hf(k) ^ [l + In (A/|fc|)] d{k) ■ a. The Berry curvature 
V X A is independent of the overall real renormalization fac¬ 
tor since it measures the “complex phase” of the Hamiltonian 
eigenstates as they are parallel transported in the BZ. The 
Chern flux through a small sphere enclosing a double-Weyl 
point remains unaltered and so do the associated topological 
quantities. 

De^ite the similarities between the present work and 
Ref.[^ the conclusions are in sharp difference. The RG fixed 
point in Ref.[35|is defined, in our convention, as the composite 
parameter ~ flowing to a flxed value with 0 and 

y/p —?► oo. In stark contrast to our results, the in Ref. [ssl 
vanishes exponentially under RG flow, which leads to the fact 
that the physical properties, such as C, k, and etc., are the 
same to noninteracting ones. 

In the end, we briefly discuss the effects of the short-range 
interactions and disorders. The Lagrangian density of a short- 
range interaction can be written similarly as 
where Ej = l 2 ,cri. A short-range coupling at tree-level is 
stronly irrelevant and scales as gj = b~^gj^R. A term 
may be generated under RG that drives the short-range 
couplings to strong coupling, similar to the situations in the 
parabolic semimetal a*^'*^ . However, due to the fact that the 
vanishes near the flxed point, the short-range couplings 
remain irrelevant and negligible. The effects of the disorders 
are more intriguing and detrimental. The Lagrangian density 
of a disorder can be written as Vjtp'^Mjijj, with Mj = I 2 , Uj. 
If we choose Gaussian white noise distribution for the 
disorder according to {{Vi{x)Vj{x'))) = lXij5^^\x — x'), 
we perform the average over disorder by employing replica 
method^ii^. The effective disorder terms mimic the four-fermi 
interactions but nonlocal in imaginary time, i.e. So ^ 
f d^xdTdT'Aj-ipKx, T)Mjt/ja(x, (f, r'), 

where a, b are replica indices. The disorder average vertices 
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are marginal, Aj = Aj^n, at the tree-level RG analysis. 
Hence, a more thorough treatment including one-loop correc¬ 
tions is needed, which we leave for the future studies. 

Note added-Dunng the journal review process, we found 
a preprint^ working on similar topic. Ref. also finds a 
new fixed point within one-loop weak-coupling RG analy¬ 
sis in the presence of the long-range Coulomb interaction, 
where the Coulomb interaction gets screened anisotropically 
and specific heat, C, receives a logarithmic correction, sim¬ 
ilar to the conclusion of the present paper. However, the 
qualitative differences between the results here and those 
in Ref. [3^ originate from the different way of introducing 
anisotropy, the second term in Eq. ([T]) involving the bosonic 
field (f) and the anisotropy variable p. Due to the subtle dif¬ 
ference of characterizing the anisotropy, the fixed points are 
different. Unlike the RG fixed point in this paper defined as 
A ~ constant, the RG fixed point in Ref. is 

defined as ~ e ^/77 —>■ constant, which leads them to the 
result that the logarithmic correction to the specific heat is 
SC ~ —T^a ln(To/T), in constrast to our result in Eq. (fTsT l. 
5C ~ ln(ro/r). Very recently, the double-Weyl 

semimetal is also proposed to be realized in SrSi 2 ^. 
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Appendix A: One-loop RG corrections for donble-Weyl 
Semimetal in clean limit 


The action of a double-Weyl fermion coupled to the long- 
range Coulomb interaction is given in Eq. ([T]i in the main texts. 
We can define the fermion Green’s function and the (boson) 
scalar Green’s function as 


Go(fc,w) 

Do{k,uj) 


iui + d ■ (7 

Lip- -b + v^kl ’ 

k'i + ilkl ’ 


(Al) 

(A2) 


where we define fc_L = (fcx, ky) and |fcj_| = fcj_ = + 

The Eig. |3a) illustrates the Coulomb interaction induced 
fermion self-energy (fc, w) 


T.-%k,u:) = -P j [ Go{qM)Do{k-q,OJ-n) 

Jn jq 

= f _(A3) 

2 Jq +^2^2 1 ^^ _ 

where we introduce the abbreviations = fj^dfl/(27r) 
and = f' d?q/ (27r)^, and the prime means the momentum 



(a) (b) 


FIG. 2. Feynman diagrams for the self-energy corrections due to 
the long-range Coulomb interaction. The red curvy lines are boson 
propagators introduced for performing Hubbard-Stratonovich trans¬ 
formation of the four-fermion Coulomb interaction. The blue lines 
are fermion propagators. The boson-fermion vertex is —iel 2 - 


integral within a momentum shell between [Ae“^,A], with 
£ <C 1. We adopt the RG scheme introduced by B.-J. Yang 
et al. and introduce the large momentum cut-off A along q± 
direction, while there is no restriction for the integral along q^. 
We take the calculation for the correction to d^a^ for exam¬ 
ple. We consider Tr[(T^i9fe^E“(fc, 0)]/rr[cr^cr^]|jr^p, which 
gives the correction to d 3 {k). 


Tr[a^dk^E^^{k, 0 )] 



VveH 1 

27r2mu, X ^ V 22 + A2 (^2 + 1)^ 


^ ^e2 
47r2 


(A4) 


where we introduced z = ^ryg^/A, and dimensionless A = 
y/rjA/{mvz). If we solve the RG equations, we will see that 
the parameter A is irrelevant and flows toward zero, and we 
show the leading contribution in the last line above. 

For the correction to d 2 cr'^ for example. We consider 
Tr[aydkjky^^^{k,0)]/Tr[ayay]\^^^, which will give the 
correction to ^2 (fc). 


Tr[aydkjk,^^^{k,0)] 

Tr[ayay] 

^ 8076 ^ f' __ 

m Jq ^m-‘^q]_ + vlql {q]_ + -qqlf 

=^ t r _ 

47r ni JAe-^ J-oc y^m“2g4 _|_ ^ 2^2 (-^2 _|_ J^g, 2 ^ 


3 


2 TT‘^mVz 




A^i 


2Tr'^mVz Jo ^ + ^2 (z2 + 1)3 

Cy/rje'^i 
2 'K'^mvz ’ 


(A5) 
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where c = ln(3 + 2-\/2)/2 ~ 0.881. During the calculation, 
we introduced the momentum cut-off for since if there is no 
restriction, we will get an artificial logarithm of A. It is also 
physically intuitive to introduce a momentum cut off for q^. 
Since if we choose Ax to be the largest momentum scale and 
perform the momentum shell integral along q±, the largest 
momentum Az along qz should satisfy VzAz ~ 
which are the energies along q± and qz- Therefore, we can get 
an identity as ravz/{A\/Az) — 1. For Ax — A^, we can get 
mvz/A ~ 1 and thus A = ^A/{mVz) — From rotation 
symmetry, we know the correction to di (fc) is the same to that 
of d 2 {k). Combining the corrections with the bare terms, we 
get 


d±{k) ■ a± 


1 + 




2 1 


d 3 ,{k)c 


y/rje^P 


1 + 


ItT^Uz 


(A6) 


where we define a± = (cri,cr 2 ). It may seem that there is 
some inconsistency in performing the RG calculations. The 
more valid way to perform the calculation should be stated 
as follows. Since we set Ax to be the largest momentum 
scale and we perform momentum shell integration within 
kj-l C [Axe“^,Ax], we should consistently introduce the 
large momentum cutoff for the integrations of qz in the self¬ 
energy calculations. However, we note that these will simply 
complicate the coefficients of the corrections and the structure 


of the RG equations will remain the same, i.e. the fixed point 
structure will remain the same. 

The Fig.|2b) represents the boson self-energy n(fc, to) 


n(fc, w) = 


Tr 


Go{q, Q,)Go(q + k,fl + ut) 


(A7) 


Since the boson propagator is frequency dependent, we can 
focus on static n(fc, w = 0). After frequency integral, we get 


n(fc,o) 



1 


d{q) ■ d{q + k) 
Eq+kEq{Eq+k + Eq) 

(A8) 


where we define = m '^q'^ + v'^q'^. After expansion to 
quadratic order in k, the integrals give 


Pi 

Qtt'^Vz 


{kl + ky) + 


nPvze^i 2 
247r2A2 


(A9) 


Combining the corrections and the bare terms, we obtain 
1 


207 


[kl + ky) 


\/vP 

1 + ^ 1 1 

+ AtA 

"ilt^Vz 

2 ^ [ 


1 + 


9 9 

m VzC 
12 Py/TjA‘^ ■ 
(AlO) 


The renormalized action after inclusion of the self-energy 
corrections due to the long-range Coulomb interaction is 


Sl =J 




207 


1 + 


y/rjp 

3tt^Vz 


11 + ^03^3 


AtT^Vz 


ij}- 






I- 

We rescale the parameters as r = x = xrP, y = The RG equations are 

yne^, z = e = and tj) = 

_ 2 /2 

Z^ (f)R to bring the action back to the original form. We 
obtain 




1 , 

4Pvz 


Cy/rie^ 


APvz 


Z^ = 


Zy = 




1 + 


076^ 

3tt^Vz 


1 + 


2 2 
m VzC 

127r20yA2 

0?e^ 1 ^ 
1 -I- ——I 
3tt^Vz 


127r2ypA2 _ 

2-1-1 


(A12) 

(A13) 

(A14) 

,(A15) 


1-b 


0?e 

37r2ti2 


1 + 


9 2 

mvze 

127r207A2'' 


d In Vz 07e^f? 

= z - Zi + 


dl 
dlnm“^ 


47r2uz 


dl 

dlnp 

dl 

din e2 

dl 


= {z-2 


— 2(1 — Zi) + 


APVz , 

nPvze^ 0?e^ 


= 2 - 1 - 


127r2y^A2 Stt^Vz 
y/rje^ rPvzG^ 
Q'k'^Vz 2AiT‘^yfyA^ 


(Alb) Introducing the dimensionless parameters, 
(A17) 


a = 


1271-^A’ 


\ — 

^ — 487r^0A^ 1 


(A18) 

(A19) 

(A20) 

(A21) 


(A22) 
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we obtain the RG equations, 
dlnu 


M 

dlnm 


M 
din a 
M 

din A 
M 


If we hold Vz and m ^ fixed, we get 


3 


"' + iA' 

(A23) 

3ca2 


4 A ’ 

(A24) 

1 


= z-l-2A-- —, 

(A25) 

= z Zi — 2 — 4A. 

(A26) 


2 = 2 - 


3c or 

TT’ 


3 


3(1 — c) o? 
4 T 


(A27) 

(A28) 


The RG equations for double-Weyl semimetals in the presence 
of long-range Coulomb interaction are 


da 

~di 

dX 

Yi 


= 2A 


1 - 2A- 
l-2A- 


3c -f 2 a 


21 


4 A 
6c — 3 ct^'' 


8 A 


(A29) 

(A30) 


{koyh) and (q^)^ = (qi ± fci/ 2 )^ -|- (q 2 ± ^ 2 / 2 )^ and the 
repeated subscript indices j means summation over j = 1 , 2 . 
We also introduce the rescaled dispersion + k^. 

We then introduce q = q -t- xk so that q(q -I- k) = (q — 
xk) [q -f (1 — a;)k] and perform the integration of qg and q^. 
We obtain the static polarization function as 

n(0, \/mkj_, —) 

Vz 

[ -dj{ql)dj{ql) + x{l - x)kl 
Stt^Vz Jo jq^ a;(l - x)kl + x{q'\_Y -b (1 - x){q'^Y 

(B2) 

Now we can examine the leading terms in ko and k±. First if 
we set k± = 0 , the result after regularization is 

n(0,0, —) =Ifal oc [fcal, (B3) 

which is linear in ^ 3 . If we set k^ = 0, the integral can not 
be performed analytically. But we can factorize out the k± 
dependence to see how the result scales with k±. We find that 
the result is 

n(0,V™fc_L,0) = - fci / f{x,y) cx fci, (B4) 


We can see that fixed points are located at (a, A) = (0, 0), 
and (0, 1/2). Linearizing around these two fixed points, we 
find that the fixed point (0, 0) is the unstable Gaussian fixed 
point, and (0, 1 /2) is the stable fixed point controlled parame¬ 
ter defined as the ratio of long-range Coulomb interaction and 
the anisotropic parameter. The RG flow diagram is shown in 
Fig. 2 in the main texts. 


Appendix B: Random Phase Approximation analysis of 
screened Coulomb interaction in double-Weyl semimetals 

We use RPA analysis to examine the screened Coulomb in¬ 
teraction in double-Weyl semimetals. We will focus on po¬ 
larization function n(A:, w) illustrated in Fig. 2(b) in the main 
texts and perform the integral without restricting integrating 
range. For clarity, we relabel the frequency and the momenta 
{uj,kx,ky,kz) —?► (fco, fci, fc 2 , ^ 3 )- The polarization function 
after proper scaling of the variable is 

--^n(fco, —) 

Zme^ Vz 


where 


f{x,y) = 


4a:^ + Axy'^ + x 

(2x — 1)^ -b 4y^ 


X In 


( 21 + 1 ) 2 + 47/2 


.(B5) 


We can see from Eq. (IB4b that the leading term in k± is still 
quadratic. We can conclude that the leading terms in RPA 
analysis is 


n(0, s/mk±, 



k± + 1^21, 


(B6) 


which is consistent with the RG analysis. 


Appendix C: Large-N/ analysis 

In this appendix, we will illustrate how the infrared diver¬ 
gency near the stable fixed point arises via a simplified large 
Nf analysis. The strong-coupling analysis starts with the ac¬ 
tion at the stable fixed point. 


qo{qo + ko) - 93(93 + ko) - dj{q- ^)dj{q + |) 


[9q + i?2j _|_ /jp) 




S,= {Ho-i4>)^XNf J (9i +1931) I'CCl) 


90(90 + ko) - 93(93 + ko) - dj{q- ^)dj{q + |) 


/ q,x 


(q + xk) + x(l - x)k2 + x(9)[)4 + (1- 


x){q± 


-A 


2 ’ 


where Nf different copies of fermions are introduced and we 
suppress dimensionful numbers and the boson propagator is 
from the RPA calculation. Here instead of evaluating the 
Nf ^ 00 limit exactly, we use the RPA results that capture 
(Bl) the correct momentum dependence in each direction. Since 
^ we are only interested in how the infrared divergence appears, 

'^here ^ dx and we introduced the Feynman pa- the use of the simplifled RPA result can be justified. It is im- 

rameter x, which leads to two vectors q = ( 90 , 93 ), k = portant that the electric charge does not appear in the action 
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since it is always possible to absorb the constant into the bo¬ 
son field by redefining the field. 

Given the approximate action, we can evaluate l/Nf cor¬ 
rection. The electron self-energy with the momentum cutoff 
A in the quadratic direction is 


n \ ^ f daif^ + QW 


1 


Ek 


(C2) 


where dacr“ 


q ^k+q 

{k\ — -I- 2kik2<7^ -b ^3(7^ and the dis¬ 


persion E{k) = \J'Ylia^daY■ The correction can be read off 
by considering the fc —>■ 0 limit as 


d^f{k + q) 


dda 


k —^0 


1 






(C3) 


We can see there is a logarithmic divergence at the infrared 
limit (/i —0). Hence, the anisotropic screening induces 
the logarithmic correction similar to conventional graphene 
physics. 


Appendix D: Dynamic conductivity at noninteracting limit 

Within linear response theory, we first start from the Mat- 
subara formalism to calculate the current-current correlation 
function 


n 




fp 

= -e^ / (T, [T^r) J^(0)]), 

Jo 


(Dl) 


and perform analytic continuation to the real frequency as 

=Ilai3iiUJn ^ U!+ iri~^). (D2) 

In the end, we can extract the dynamic conductivity by ex¬ 
tracting the imaginary part 


j fill 


(w) = - 


Im [n^^(w)] 


(jJ 


(D3) 


Before the indulging in the calculations, we first note that 
due to the rotation symmetry, ^ 

^zzi'iojn)- We will discuss separately 

and Hzziiuin), which lead to the conductivity axx(^) = 

ayy(uj) z/z azz{oS). 

The current components are 

Jx = j ^ 

Jy=['^l f”—cr^ + —(tA il}q= f Jy V^„(D5) 

J q \ Vfl Tfl J Jq 

Jz= ! = [ 'fplJz fpq- (D6) 

J a J a 


fq Jq 

and the diagonal current-current correlation function in the 
Matsubara domain can be expressed as 

o2 


n^/x(*Wn) = ^ X/y [Jy.Q{q,iPm)JyQ{q,iPm + fw„)] 


where Q{q,iujn) is the noninteracting fermion green’s func¬ 
tion in the Matsubara domain as 


c;(k,ia;„) = - 


1 


iuJn + p - 710' 


(D8) 


where the 'Ho is the Hamiltonian density of the system. After 
straightforward derivation, below we list the main results. 

(A) Drude weight at zero frequency, (w) 

a;—>-0 


^xxjyy (^) 
^2 


cj—>-0 


Sttuiv- 


■(2mT) Li 2 {-e^^'^) + Li 2 {-e : 

(D9) 


O’zz(w) 

mvze 


uj—^O 

2 r 


16 


,5 . (DIO) 


The Drude weights show different behaviors at different lim¬ 
its. 

( 1 ) p/T < 1; 




C(2)e 


cj—)-0 


Sttuiv. 


,{uj) 




mvze^ 


uj—¥0 


-(2mr)5(^); (Dll) 
■ln(2)<5(^). 


(D12) 


(2) p/T » 1: 

XXjyy (^) 


a. 


at finite 


(B) Dynamic conductivity at finite frequency, tr, 

0 ): 

^XX jyy (‘^) 

^^■>0 


'-) ;(Di3) 

(D14) 


txJ>0 

—- (muj) 

24:TrmVz 


tanh(^^^^) +tanh(^^^^) ; 

(D15) 


Ckzzioj) 

u ;>0 

= L„h ('±^'\ tanh ()1 

128 [ \ 4T ) V 4T ) 

Combining both the zero frequency and finite freqr 
the dynamic conductivity can be expressed as ax 
(D7) cryy(w, T) = 3^^(2mT)$ I (w/T. u/T), anH ^ 


(D16) 

requency and finite frequency parts, 
can be expressed as axx{oJ,T) = 
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<i>a(a;/T, /^/T) with the scaling function 


$_L(a, b) = - [Li 2 {-e^) + Lz 2 (-e '*)] 6 {a) + 

1 


16 


^ ,a b. , .a b. 

tanh(- + -) +tanh(- - -) 


^z{a, b) = - [Lii{-e^) + Lii{-e '’)] S{a) + 


1 

--Cl 


, ,a b. , ,a b, 

tanh(- + -) +tanh(- - -) 


5(D17) 


. (D18) 


Appendix E: Diamagnetic susceptibility at noninteracting limit 

In the noninteracting limit, the Hamiltonian is 

~2 „2 c\ 

jj- Q-x^y ^QxQy /triA 

Hows = - CTx H- CTy + v^qziJz- (El) 

m m 

In order to calculate the diamagnetic susceptibility, we will 
use the Fukuyama formula as 

Xd = X] / [QolaQolbQolaQolb ], (E2) 

where Qq is the fermion Green’s function in the Matsubara 
domain, the n summation represents the Matsubara frequency 

I 


sum, and 7 ^ = dH/dqa, with a being the direction axis that 
perpendicular to the direction of the magnetic held. Before we 
go into the calculations, we hrst note that due to the anisotropy 
it is expected that the diamagnetic susceptibilities for the cases 
with B = Bx and B = Bx are the same. However, the dia¬ 
magnetic susceptibility for the case with B = Bz should be 
different to the two former cases. Let us discuss each case 
separately below to see the temperature dependence of the dia¬ 
magnetic susceptibilities in different cases. Below, we set the 
chemical potential to be zero. 

First, we choose B = Bx and the result should be the same 
to the case of B = By. Now, we have 7 ^ = —{2qy/m)ax + 
{2qx/may) and 7 ^ = The Fukuyama formula gives 


X / [^o7y^^o7^^^o7y5o7^] (E3) 


1 

m? (3 


Y. h’- 


q^Qocy yQ^a zQoCTyQocj z + 


-f-qyQo<^xQo<^ zQo'y xQo<^ z 


±,(/) , _L.(7/) 

= Xd +Xd ■ 


(E4) 


After expansion and exchanging q^ and qy for Xd 

= Xd’^^^- After performing the Matsubara frequency 
summation, we get 


Xd = 


16e^v^ 


m 


tanh(^) sech'^{-^y 


Sq^d^d^ 


sech?{i^) tanh^(^) 


192d4r3 


sec/i^(^) tanh(^) ^ bsechr[^) 


5tanh(^) 


32d5T2 


64d6r 


32d7 


■ (E5) 


The momentum integral is complicated, but since we are only After the rescaling and straightforward algebra, we get the 
interested in the temperature dependence, we can factorize out diamagnetic susceptibility for B = Bx, 
the temperature dependence by rescaling 

qx = {mT)^x, qy = {mT)iy, qz = Tv~'^z. (E 6 ) 

I 


Xd — 


-e^Vz 


dxdydz 2 

’ 2 tanh(|) ^ sec/i^(|)’ 

/-oo (27r)3 [ 

[ d3 d 2 J 


-f 


+ -X^^2^3 

o 


sech?{^) tanh^d) 


-f 


-e^Vz 


Gsech^ (I)tanh(|) 

52 


15sech^(|) 


30tanh(|) 


d7 


(E7) 


where we introduce d = (df -f d^ + d§) , and di = — 

y^, d 2 = 2xy, and da = 2 . The result above should be 
the same for the case with B = By. Now let us check 
the temperature dependence for the diamagnetic susceptibil¬ 


ity in the presence of 73 = Bz. In this case, we need 

lx = i^qxl'm)ax + {2qyjm)ay, and Xy = -i2.qylm)(7x + 
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{2qx/m)ay. The Fukuyama formula gives 


Xd 




QoXxGoXy^oXxQoXy 


(E8) 


32e^ 


Xd =- 


Iq L 


g^tanh(^) qj_sech'^ 


4d3 


8 d‘^T 


After expansion and performing Matsubara frequency sum¬ 
mation, we get 


256e2 

m® 




sech?{i^) tanh^(^) 


sech?‘{i^) tanh(^) 
32d5r2 


5sech^{^) 5tanh(^) 
64d6r 32d7 


(E9) 


We can again factorize out the temperature dependence. We find 

I 


Xd = 


4e"i 


axaya^ 


\ z ' 


mvz (27r)3 [ rfs 

4e^r dxdydz , 2 


J2 


ri+4(^2-2/2)4 


sech'^ii) 6sec/i2(^) tanh(l) 

-=-^ tanh2(-) H-1^4-^ + 

d^ ^ 2 ^ dP 

lbsech?{^) 30tanh(|)' 


mvz 


(ElO) 


Hence, in the presence of i? = Bz the diamagnetic suscep¬ 
tibility is actually linearly proportional to the temperature. 
Combining the results of the cases of i? = Bt''± and B = Bz, 
we expect that the diamagnetic susceptibility in the presence 
of magnetic field in arbitrary direction B = Bf should show 
the temperature dependence as 

Xd ^ (sin2 9x0 + cos^ 0 T) , (Ell) 


where xo is a constant independent of T and we introduce the 
periodic function with 9 being the angle between the magnetic 
field and the z—axis, i.e. B-z = B cos 9. The square of the pe¬ 
riodic functions roughly gives the correct angular dependence 
with xd( 6 ) = xd{9 + 'x). Therefore, at low temperature limit 
T —0 we expect that the constant diamagnetic susceptibility 
dominates. 
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